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1 Introduction 


Beginning from the first identification of extragalactic source 3C 273 sd y-ray astronomy 
has made great progress due to unprecedented detection capability of the new generation of 
telescopes” ”. The large area telescope (LAT) onboard Fermi satellite has detected more 
than 3 000 y-ray sources after first four years of operation, out of which about half are 
extragalactic sources”. The extragalactic y-ray sky is dominated by emission from blazars 
(about 98%), which is an extreme subclass of active galactic nuclei (AGN) with relativistic 
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jets pointing close to our line of sight. Blazars have been classified as BL Lac objects 
(BL Lacs) or FSRQs according to whether the rest-frame equivalent width of their broad 
emission lines was greater than (for FSRQs) or smaller than (for BL Lacs) 5 A. Due to very 
high luminosity and violent variability, their y-ray emission is believed to be produced in 
the relativistic jets making a small angle to our line of sight. The intensity is significantly 
boosted by the Doppler beaming effect because of light aberration through transformations 
I, = I, with 6 = v/v’ from intrinsic to observational frames, with the primed quantities 
referring to values measured in jet co-moving frame, while 6 = 1/I’(1—{ cos 0) is the Doppler 
beaming factor. 

The study of y-ray emission of blazars plays an important role in understanding the 
jet physics (jet launching, energy dissipation and particle acceleration), the cosmic evolution 


B, & B 


history and the origin of ultra high energy cosmic rays . The y-rays can be explained as 


inverse Compton (IC) emission of accelerated electrons in the jet. The IC emission in BL Lacs 
arises due to internal synchrotron photons and called synchrotron self-Compton (SSC)™. 
The multi-wavelength catalog of blazars, Roma-BZCAT® v5.0, is the most comprehensive 
catalog in the literature, which contains 3 561 confirmed blazars or candidates”. Recently, 
Fermi Fourth LAT AGN Catalog (4LAC) presented a sample of nearly 3 000 blazars or blazar 
candidates at high galactic latitude’. Therefore, there may be >500 blazars having not been 
detected with y-ray emission (y-ray quiet), although they have similar emission properties at 
lower energy wavelength bands (radio through X-ray) compared to y-ray detected blazars. 
The origin of this discrepancy is not yet known. Many studies have focused on this question 
and found that y-ray loud blazars have relatively large jet opening angle, larger apparent 
superluminal velocity and higher brightness temperature relative to y-ray quiet blazars “7 
All these properties are attributed to a relativistic geometric effect - Doppler beaming effect . 
The smaller viewing angle in y-ray loud blazars makes apparent larger jet opening angle and 
larger Doppler beaming factor, which give rise to larger apparent superluminal velocity and 
higher brightness temperature. In addition to this geometric effect, it is expected that the 
origin of the discrepancy between y-ray loud versus y-ray quiet blazars might be intrinsic. 
However, it is not yet known. 

The spectral energy distribution (SED) of synchrotron emission from blazars is usually 
curved and presents a significant bump in v — vf, frame, which peaks between infrared to 
X-ray bands”. The curvature of broadband SED of the synchrotron bump is different 
among various blazars, which means that some SED bumps are broad while others are 
narrow. If the SED (broad band or single band) is fitted by a log-parabolic function, 
lgv fy = lg vp fu, —b(igv—lg vp)’, the coefficient of the second order term b can be considered 
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a “surrogate” to measure the curvature. The curvature is found to be negatively correlated 
with the SED peak frequency, which offers supporting evidence of a statistical or stochastic 
particle acceleration at work in blazar jets, through studying the broad band observations 
or single band (X-ray) observations. It should be noted that the broadband curvature of 
a blazar is not changed with the transformation from observer to AGN frame (Cosmological 
redshift) and jet co-moving frame (Doppler beaming). However the single band curvature 
is effected due to Doppler effect causing the spectral shifts. In this work, we investigate the 
possible relation between y-ray emission and the SED curvature for BL Lacs. The curvature 
in SED is produced due to an intrinsic curvature in the emitting electron energy distribution 
(EED) arising due to competition of particle acceleration and cooling. The Section 
describes our sample and Section B presents the methods and a discussion on our results. 
We summarize our findings in Section M. We assume a ACDM cosmology with values from 
the Planck results in our calculation; in particular, Qm = 0.32, Qa = 0.68, and the Hubble 
constant Hy = 67 km-s~!-Mpe7!™™ 


2 The Sample 


Starting from Mrk 421”, Massaro et al.” found that curvature is an important feature 
of BL Lacs in X-rays. Even the Fermi-LAT y-ray spectra of bright blazars show curvature. 


The high energy photon spectra can be fitted by a log-parabolic law a 
E —a—ßlg(E/Eo) 
F(E)= K| = ; (1) 
Eo 


with spectral index a at reference energy Eo and £ is the curvature of bump. The corre- 


sponding peak energy of the spectral bump becomes 
E, = Bylo") . (2) 


Ackermann et al.” found that the LAT spectra of about 160 blazars show deviation from a 
power-law F(E) = kK E~” and fitted with log-parabolic law. This suggests that the curvature 
is a common feature of blazar SED in any energy band. 

The Fermi-LAT detected 1591 blazars at high galactic latitude in its Third LAT AGN 
Catalog (3LAC) based on first 4 years of operations”. Based on 3LAC, Fan et al.” collected 
available spectral data of 1392 blazars from NEDY to build broadband SEDs of synchrotron 
bump and fitted them by a log-parabolic function. We compile the synchrotron spectral 
curvature (b), synchrotron peak flux (F, = vy) f,,), and integrated y-ray flux F, in 0.1~100 
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GeV band for confirmed 620 BL Lacs in 3LAC™. Wu et al. presented the Doppler factors 
of a BL Lacs. Cross matching two catalogs, we compile the spectral parameters and Doppler 
factors for a sample of total 170 BL Lacs. 


3 Methods and Results 


A curved SED necessary suggests that the underlying emitting EED may be curved. 
A log-parabolic EED can arise due to stochastic nature of acceleration gain. Anjum et 
al." found that only BL Lacs show the signature of stochastic acceleration, whereas FSRQs 
do not show any signature of such acceleration mechanism. The broadband SED in the 


log-parabolic model can be described as 


lev fe = lg vfu — blgv — lg ws) (3) 


where b measures the width (curvature) of the SED bump. The synchrotron bump is par- 
ticularly important and reveals the intrinsic curvature of EED. Since the accretion disk 
continuum emission in FSRQs are usually presented as a big blue bump at ultraviolet wave- 
length band, affecting accurate measurement of the curvature parameter, we only selected 
BL Lacs for purpose of our study. The selection of BL Lacs is further motivated by the fact 
that the high energy bump in FSRQs arises due to external Compton (EC) emission, there- 
fore, their curvature of IC bump may be related to the complex seed photon distribution 
and any absorption”. Although the curvature of high energy bump in BL Lacs mimics the 
synchrotron curvature as it arises due to internal synchrotron seed photons”, it is hard to 
constrain the intrinsic of the EED from the IC bump, since the observed SED curvature of IC 
bump may also depend on the IC cooling regime. Paggi et al." found that the IC curvature 
in Thomson regime is smaller than in Klein-Nishina regime. Therefore, we constrain the 
intrinsic curvature of EED r from synchrotron spectral curvature b rather y-ray curvature 
B. The intrinsic curvature of EED r ~ 5b and y-ray emission in BL Lacs is expected to be 
intrinsically related and can be used to investigate the origin of discrepancy between Fermi 
detected BL Lacs (FBLs) and non-detected ones (NFBLs). 
3.1 FBLs versus NFBLs 

A large amount of multi-wavelength data for the objects in the Metsahovi radio ob- 
servatory BL Lac sample was collected by Nieppola et al, which is supposed to have no 
selection criteria (other than declination) in addition to the ones in the original surveys. 
The Metsahovi radio observatory BL Lac sample includes 381 objects selected from the 
Veron-Cetty & Veron BL Lac Catalog™, and 17 objects from the literature, of which many 
sources are from the well-known BL Lac samples like 1Jy, $4, S5, Einstein Medium Sensi- 
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tivity Survey (EMSS), Einstein Slew Survey, and ROSAT Deep X-ray Radio Blazar Survey 
(DXRBS). Based on the multi-wavelength data, the SED of each source were constructed 
in the lgv — lgvF, representation, of which the synchrotron bump was fitted with a log- 
parabolic function. Based on this sample, Wu et al." collected the available data at 330 
MHz, 360 MHz, 408 MHz, and 1.4 GHz from the Astrophysical Catalogues Support System 
(CATs) maintained by the Special Astrophysical Observatory, Russia, and also the available 
VLA or MERLIN core and extended flux, resulting in a sample of 170 BL Lacs. The low fre- 
quency radio power can be a reliable indicator of the intrinsic radio power, while the Doppler 
beaming can affect the observed radio power of the core Therefore, with the VLA or 
MERLIN core and the 408 MHz luminosity and assuming a jet speed I’ = 5, consistent with 
Metsahovi radio monitoring studies ™, the Doppler beaming factors have been estimated ™. 
We compiled the Doppler factors and curvature parameters of these 170 BL Lacs. By cross 
correlating this sample with 3 LAC”, 34 out of these 170 BL Lacs are not found to have 
y-ray emission (NFBL). These NFBLs are: NRAO 5, NPM1G +41.0022, 1ES 0145+138, 
MS 0158.5+0019, MS 0257.9+3429, RXS J0314.0+2445, S5 0454+84, MS 0607.9+7108, B3 
0651+428, 4C 22.21, RXS J0916.8+5238, B2 0927+35, RGB J0952+656, 1ES 1044+549, 
1ES 1212+078, 1ES 1255+244, 1ES 1320+084N, MC 1400+162, MS 1407.9+5954, RGB 
J1427+541, MS 1443.5+6349, RXS J1516.7+2918, MS 1534.2+0148, RXS J1602.2+3050, 
RXS J1644.2+4546, RGB J1652+403, B3 1746+470, RXS J1750.0+4700, RGB J1811+442, 
PKS 2254+074, Q J2319+161, 1ES 2326+174, MS 2336.5+0517 and MS 2347.4+1924. The 
remaining 136 BL Lacs are, therefore, considered FBLs. Although many of these NFBLs 
might be detected in the 4LAC, however their y-ray flux must be lower than the FBLs. 


The SED of blazars is usually curved even in a single energy band. The Doppler factor 


changes cause the peak vp = 1/6 shifts in the SED (vI, = 6*v'T/,) which affects the spectral 
curvature of blazars in Fermi-LAT y-ray band. Other than the peak shifts, the intrinsic 
difference of curvature may account for discrepancy between FBLs and NFBLs. We reduce 
the effect of the relativistic Doppler beaming and directly compare FBLs and NFBLs. We 
searched the literature for published and archival observations that allow us obtain both 
curvature parameters and Doppler beaming factors of our Fermi BL Lacs. We plot curvature 
parameter (in 1/b) against 6 for FBLs and NFBLs in Figure ET, which shows that both FBLs 
and NFBLs overlap but the curvature of FBLs would be on average smaller than that of 
NFBLs. We binned the data in lg(ô) and the find the average value of curvature in each bin, 
represented by big circles in Figure B.I. Due to only one NFBL in the last bin, we consider 
it as the average curvature value. FBLs, on average, seem to have smaller curvature (higher 
value of 1/b) and, thus, broader SED than NFBLs. 

In order to eliminate effect of Doppler beaming as much as possible, we artificially 


defined a new parameter, b — lg6, which roughly measures the curvature given the same 
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Ig(1/b) 


A 
-0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 


Notes: The dashed lines represent 4 bins of Doppler factor 6. The small points represent the values 
of individual sources, while the large points show the average value of curvature in each bin. The 


error bars correspond to one standard deviation of the average value. 


Fig. 1 The curvature and Doppler beaming effect in BL Lacs 


beaming factor. Figure BI shows the distributions of b — lgô and cumulative fraction of 
FBLs and NFBLs. It can be seen that statistically FBLs have relative smaller values of 
b — lgô than that of NFBLs, which confirm the above finding that smaller curvature of 
synchrotron SED bump is attributed to larger y-ray power. The Kolmogorov-Smirnov (KS) 
test yields the significance level probability for the null hypothesis that FBLs and NFBLs 
are drawn from the same distribution P = 5.87 x 1075, and the maximum separation of 
the two cumulative fractions is Dgs = 0.43. The TeV BL Lacs (TBLs) with synchrotron 
peak frequency vp > 10'° Hz showed a similar behaviour”. Massalo et al.” found X- 
ray band curvature of TBLs to be systematically smaller than that of non-detected at TeV 
energies. They suggested that X-ray flux can be used as a predictor for TeV y-ray detection 
of BL Lacs. However the TeV spectral curvature may not correspond to X-ray curvature, as 
the TeV emission is significantly attenuated by extraglactic background light (EBL). Since 
GeV 4-ray emission is not absorbed, we argue that it can provide better constraint on TeV 
detection as compared to X-rays. 
3.2 SED Curvature and y-ray Dominance 

The fact that synchrotron SED is significantly curved inevitably implies a curvature in 
the steady electron energy distribution N(y). This curvature may be a result of particle 
cooling and acceleration which might be relevant on different times. Based on scenario of 


re-acceleration rather than continuous injection, the curvature can be related to stochastic 


88 RM +E 39 & 


Notes: The Kolmogorov-Smirnov (KS) test yields the significance level probability for the null 
hypothesis that FBLs and NFBLs are drawn from the same distribution P = 5.87 x 1075 and the 


statistic (the maximum separation of the two cumulative fractions) Dgs = 0.43. 


Fig. 2 The distribution and the cumulative fractions of FBLs versus NFBLs 


acceleration term in the Fokker-Planck kinetic equation". A more efficient acceleration 
of stochastic type makes energy distribution electrons relatively broader, manifested by a 
smaller curvature ™. A higher acceleration efficiency implies a lesser time (t = R/c) spent in 
acceleration region by the particles and large number of acceleration steps, which corresponds 
to a smaller size of the acceleration/emission region”. Therefore, the spectral curvature 
of EED is expected to be proportional to the size of blazar region. As the y-rays from BL 
Lacs are believed to be SSC emission of the same non-thermal electron population emitting 
low energy synchrotron emission, a smaller size would enhance the synchrotron seed photons 
energy density us = L,/47tR?c6* in the jet, that would lead to higher Compton dominance 
CD = u,/up, i.e, a higher relative y-ray power. This implies that among the BL Lacs with 
the same synchrotron luminosity Ls, the source with relatively high y-ray dominance may 
have relatively compact emission region. 

Abdo et al” provided broadband simultaneous or quasi-simultaneous spectral data 
from radio through y-rays, of 48 LAT Bright AGNs, within the first 3 months operation 
of Fermi-LAT, from 2008 August 4 to October 31. The broadband spectral data are de- 
rived from many ground-based and space-based observatories, including Swift (UVOT, XRT, 
BAT), Effelsberg, OVRO, RATAN, GASPWEBT, Spitzer-MIPS and AGILE. We show an 
example of SEDs of two BL Lacs (PKS 0048-09 and S5 0716+714), which are presented 
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in Figure BA. These sources have similar synchrotron luminosity and peak frequency but 
different y-ray power, i.e., the Compton dominance. Their quasi-simultaneous broadband 
SEDs are compiled from reference[42] shown as large and small black squares and magenta 
circles for PKS 0048-09 and S5 0716+714, respectively. The S5 0716+714 included 3 obser- 
vations at optical and X-ray bands during those three months. We use the average value 
of fluxes for this BL Lac object. 


Ig(vL,/10J-s7) 


= PKS 0048-09 
e $5 0716+714 


8 10 12 14 16 18 2 22 24 26 28 
lg (v/Hz) 


Notes: Quasi-simultaneous data points are presented as large squares (PKS 0048-09) and circles 
(S5 0716+714), while the small ones show the non-simultaneous data. The solid lines represent 
the synchrotron+SSC model curves. Both sources have same jet parameters other than size R and 


curvature r assuming r x R. 


Fig. 3 Broadband SEDs of PKS 0048-09 and S5 0716+714 


We use standard one-zone SSC model” to fit their SEDs. One-zone model is widely 
used in blazar SED modeling” I Since the coordinated variability in different wavelength 
bands is often seen (although not always) in blazars, a one-zone model assumes that bulk 
of the non-thermal emission are produced from a “one-zone” region, mostly assumed to 
be a spherical blob with radius R embedded in a homogeneous but tangled magnetic field 
B and filled with emitting electron population. Since the considered emitting region is 
always compact, its synchrotron self-absorption frequency is always large, therefore the model 
cannot account for the radio flux at observed frequencies smaller than a few hundreds GHz, 
which are produced by the superposition of several larger components”. The emitting 
relativistic blob moves with angle 0 with the line of sight, yielding a Doppler beaming 
factor 6. In order to keep accordance with log-parabolic shape of synchrotron SED bump, 


the emitting particles are assumed to be leptons with energy distribution following a log- 
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-3 
N(= Mo (=) 10-7 180/%)? (4) 
Yo 
where r measures the curvature of number distribution of electron energy". As discussed 
above, the curvature may be proportional to the size of emission region in the SSC model. 
Therefore, in our calculation, the curvature parameter (8) is assumed to be proportional to 
the radius of the emission sphere, i.e., r x R. The calculated SEDs are also presented in the 
Figure BA. It seems that both SEDs can be well fitted with same jet parameters except for the 
radii of emission sphere Rooag = 2.8 x 1018 cm and Ro716 = 5.0 x 1016 cm for PKS 0048-09 and 
S5 0716+714 (and therefore the curvature roo4g = 0.5 and ro716 = (Roz16/ Roos) rooas = 0.89) 
respectively. The other constant jet parameters are the magnetic field strength B = 0.15 G, 
the Doppler factor 6 = 24, the peak electron energy yo = 8660 and the normalized total 
electron numbers No = 6.8 x 10*°. This modeling example illustrates that a higher y- 
ray luminosity intrinsically accompanies a smaller curvature, which may be related to a 
smaller size of emission region and the processes governing the jet micro-physics, i.e., particle 
acceleration and cooling. This is further consistent with the fact that bright y-ray blazars are 
more (rapidly) variable than y-ray weak blazars from studying y-ray and optical variability 
of a sample blazars = , possibly due to smaller size of emission region and/or larger Doppler 
beaming factor in bright sources. Figure shows that among the BL Lacs with a given 
synchrotron peak luminosity and frequency, the object with the highest CD might have 
smaller curvature (i.e., hard spectrum). A high CD coupled with a lower curvature favor 
the detection of a blazar by Fermi. Therefore, we argue that CD and curvature both are 


important parameters for Fermi detection. 


4 Summary 


We investigate why are some BL Lacs detected having y-ray emission by Fermi but 
others not. We find that the SED curvature and y-ray dominance of BL Lacs might be 
intrinsically related. We select 170 BL Lacs with synchrotron SED curvature b and Doppler 
factor 6 reported in literature, and divide them into FBLs and NFBLs. We find that FBLs 
have smaller curvature than that of NFBLs even after getting rid of the beaming effect. 
We show an example that two Fermi BL Lac objects PKS 0048-09 and $5 0716+714 have 
similar synchrotron peak frequency and luminosity but different Compton dominance. The 
PKS 0048-09 shows relatively high y-ray dominance as compared to S5 0716+714. Within 
a one-zone SSC model, we find their quasi-simultaneous SEDs can be well fitted by same 


physical jet parameters except for the size of emission region and EED curvature (assuming 
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curvature being proportional to the size). The PKS 0048-09 manifest relatively smaller 
source size and a smaller curvature as compared to $5 0716+714. These results imply 
that the difference in curvature might be related to y-ray dominance and may account for 
the intrinsic discrepancy between FBLs and NFBLs. A broader SED with compact jet size 
demands an efficient stochastic acceleration. As the emitting particles have less time to spend 
in a compact jet, the acceleration gain or number of acceleration steps should be large. Thus, 
a more efficient stochastic acceleration in compact jets of FBLs makes the EED relatively 
broad, yielding a smaller value of curvature parameter. However, a smaller source size, in 
turn, implies a higher synchrotron photon energy density in the jet, producing powerful y-ray 
emission through SSC process. Therefore at a given synchrotron luminosity, the FBLs may 
have intrinsically smaller curvature and higher Compton dominance as compared to NFBLs. 
A study of large sample in future may be necessary to further investigate the relationship of 


curvature and y-ray dominance of blazars. 


SB Eh: 


Swanenburg B N, Bennett K, Bignami G F, et al. Nature, 1978, 275(5678): 298 
Madejski G G, Sikora M.ARA&A, 2016, 54: 725 

Massaro F, Thompson D J, Ferrara EC. A&A Rev, 2015, 24: 2 

Acero F, Ackermann M, Ajello M, et al. ApJS, 2015, 218(2): 23 
Ackermann M, Ajello M, Atwood W B, et al. ApJ, 2015, 810(1): 14 

Plaga R. Nature, 1995, 374(6521): 430 

MAGIC Collaboration, Albert J, Aliu E, et al. Science, 2008, 320(5884): 1752 
Kalashev O E, Kusenko A, Essey W. Phys Rev Lett, 2013, 111(4): 041103 
Abdo A A, Ackermann M, Ajello M, et al. ApJ, 2011, 736(2): 131 

Massaro E, Giommi P, Leto C, et al. A&A, 2009, 495(2): 691 

Massaro E, Maselli A, Leto C, et al. Ap&SS, 2015, 357(1): 75 

Ajello M, Angioni R, Axelsson M, et al. ApJ, 2020, 892(2): 105 

Linford J D, Taylor G B, Romani R W, et al. ApJ, 2012, 744(2): 177 
Pushkarev A B, Kovalev Y Y.A&A, 2012, 544: A34 

Piner B G, Pushkarev A B, Kovalev Y Y, et al. ApJ, 2012, 758(2): 84 

Wu Z, Jiang D, Gu M, et al. A&A, 2014, 562: A64 

Lister ML, Aller M F, Aller H D, et al. AJ, 2016, 152(1): 12 

Chen L, Bai J M. ApJ, 2011, 735(2): 108 

Ghisellini G, Tavecchio F, Maraschi L, et al. Nature, 2014, 515(7527): 376 
Ghisellini G, Tavecchio F. MNRAS, 2015, 448(2): 1060 

Chen L. ApJ, 2014, 788(2): 179 

Chen L. ApJS, 2018, 235(2): 39 

Massaro E, Perri M, Giommi P, et al. A&A, 2004, 413: 489 

Kardashev N S. Soviet Ast, 1962, 6: 317 

Tramacere A, Massaro E, Taylor A M. ApJ, 2011, 739(2): 66 

Planck Collaboration, Ade P A R, Aghanim N, et al. A&A, 2014, 571: A16 


© OND Te Wn eB 


NUN NIN ION N ON TB OE OR RR RE I e 
a ok U Ne OO DAN DoH wn Fr Oo 


92 RM HE 39 & 


27| Massaro F, Tramacere A, Cavaliere A, et al. A&A, 2008, 478(2): 395 
28] Massaro E, Perri M, Giommi P, et al. A&A, 2004, 422: 103 

29] Fan J H, Yang JH, Liu Y, et al. ApJS, 2016, 226(2): 20 

30] Wu Z, Jiang D R, Gu M, et al. A&A, 2007, 466(1): 63 

31] Anjum M S, Chen L, Gu M. 2020, ApJ, 898: 48 

32] Massaro E, Tramacere A, Perri M, et al. A&A, 2006, 448(3): 861 

33] Paggi A, Massaro F, Vittorini V, et al. A&A, 2009, 504(3): 821 

34] Nieppola E, Tornikoski M, Valtaoja E. A&A, 2006, 445(2): 441 

35] Veron-Cetty M P, Veron P. A&A, 2010, 518: 110 

36] Giovannini G, Cotton W D, Feretti L, et al. ApJ, 2001, 552(2): 508 
37] Giroletti M, Giovannini G, Taylor G B, et al. ApJ, 2004, 613(2): 752 
38] Lahteenmaki A, Valtaoja E. ApJ, 1999, 521(2): 493 

39] Massaro F, Paggi A, Elvis M, et al. ApJ, 2011, 739(2): 73 

40] Massaro F, Paggi A, Errando M, et al. ApJS, 2013, 207(1): 16 

41] Massaro F, Paggi A, Cavaliere A. ApJ, 2011, 742(2): L32 

42] Abdo A A, Ackermann M, Agudo I, et al. ApJ, 2010, 716(1): 30 

43] Yan D, Zhang L, Yuan Q, et al. ApJ, 2013, 765(2): 122 

44] Zhang J, Sun X N, Liang E W, et al. ApJ, 2014, 788(2): 104 

45] Kang S J, Chen L, Wu Q. ApJS, 2014, 215(1): 5 

46] Chen L. ApJ, 2017, 842(2): 129 

47] Blandford R D, Kénig] A. ApJ, 1979, 232: 34 

48] Konigl A. ApJ, 1981, 243: 700 

49] Vovk I, Neronov A. ApJ, 2013, 767(2): 103 

50] Hovatta T, Pavlidou V, King O G, et al. MNRAS, 2014, 439(1): 690 
51] Wu Q, Zhang B, Lei W H, et al. MNRAS, 2016, 455(1): L1 


The Curvature of Spectral Energy Distribution and y-ray 


Dominance of Fermi BL Lac Objects 
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Abstract: The extragalactic y-ray sky is dominated by blazars and their study plays an 
important role in understanding jet physics, cosmic evolution history and origin of ultra high 
energy cosmic rays. In this work, we study a large sample of BL Lac objects to investigate 
why some sources are detected in y-rays, while others not. We selected 170 BL Lac objects, 
with measured synchrotron spectral curvature and Doppler factor, and divided them into 
Fermi-LAT detected (FBLs) and non-detected (NFBLs) sources. We show that FBLs have 
smaller curvature than NFBLs, even after getting rid of Doppler beaming effect. The BL 
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Lac objects PKS 0048-09 and S5 0716+714 have similar synchrotron peak frequency and 
luminosity but different y-ray dominance and their quasi-simultaneous broadband spectral 
energy distributions (SEDs) can be well fitted by a log-parabolic synchrotron self-Compton 
(SSC) model with same jet parameters except for the curvature and source size, assuming 
curvature being proportional to the size of emission region. Our results imply that for a 
given synchrotron luminosity, the different SED curvature and Compton dominance may 
account for the discrepancy between FBLs and NFBLs. We discuss these results in context 


of stochastic particle acceleration and radiation mechanisms. 
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